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Abstract. The initial boundary- value problem (IB VP) and the Cauchy problem for 
the Kuramoto-Sivashinsky equation 

I't ~l~ ^xxxx H~ V X x — 2 ") x 

and other related 2mth-order semilinear parabolic partial differential equations in one 
dimension and in 1 N are considered. Global existence and blow-up as well as L°°-bounds 
are reviewed by using: 

(i) classic tools of interpolation theory and Galerkin methods, 

(ii) eigenfunction and nonlinear capacity methods, 

(iii) Henry's version of weighted GronwalPs inequalities, 
(vi) two types of scaling (blow-up) arguments. 

For the IBVPs, existence of global solutions is proved for both Dirichlet and "Navier" 
boundary conditions. For some related 2mth-order PDEs in R N x R+, uniform bound- 
edness of global solutions of the Cauchy problem are established. 

As another related application, the well-posed Burnett-type equations 

v* + (v ■ V)v = -Vp- (-A) m v, divv = in 1^ x R + , m > 1, 

are studied. For m = 1, these are the classic Navier-Stokes equations. As a simple 
illustration, it is shown that a uniform L p (R JV )-bound on locally sufficiently smooth 
v(x,t) for p > 2m-i i m P ues a uniform L°°(R JV )-bound, hence the solutions do not 
blow-up. For m = 1, N = 3, this gives p > 3 reflecting the famous Leray-Prodi-Serrin- 
Ladyzhenskaya regularity results (L p ' q criteria), and re-derives Kato's class of unique 
mild solutions in M. N . Bounded classic L 2 -solutions are shown to exist for TV < 2(2m — 1). 



1. Introduction and motivation 
The role of the Kuramoto-Sivashinsky equation 
(1.1) v t + v xxxx + v xx = ~ (v 2 ) x 
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is well known in contemporary nonlinear mechanics and physics. 

This equation arises as a model in hydrodynamics (a thin film flow down an inclined 
plane in the presence of an electric field), in combustion theory (propagation of flame 
fronts), phase turbulence and plasma physics, as well as one model for spatio-temporal 
chaos and in many other physical phenomena; see [71] for a nice short review of applica- 
tions with key original references. 

First results on global existence of classical solutions of (11.11) go back to the 1970s (cf. 



E5]) and the 1980s; we refer to El QH QU ESI US S31 S3 S3 EH EH E21 [H EQ] , where 




further important references can be found. A large part of previous studie of (11.11) was 
devoted to periodic problems. See HJJ for important references from the 1980s and 
also to some classes of particular solutions and the local behavior of solutions, where very 
interesting results (for instance, chaotic behaviour, estimates of dimension and structure 
of attractors, bifurcation theorems, etc.) have been obtained. Other papers with deep 
results are devoted to existence of periodic solutions and traveling wave solutions with 
complicated dynamics. Numerous contributions are dealing with dynamical analysis of 
the one dimensional Kuramoto-Sivashinsky equation. 

The present paper is devoted to a review of approaches that lead to global existence 
and blow-up for the Kuramoto-Sivashinsky equation and modified versions of it in one 
dimension and on M. N , the latter being much less developed in the literature. 

One of the key mathematical features of the KS-type PDEs is that an a priori L 2 -bound 
of solutions v(x,t) of the form 



is straightforward, while the main question is how to use this information to get stronger 
estimates in Sobolev spaces and eventually in L°°. 

To this end, we shall use some known methods and compare their strength with the 
new techniques to be reviewed and developed throughout this paper: 

(I) In Section[2]we employ classical interpolation and Galerkin methods for study global 
existence of the initial-boundary value problems (IBVP) in one dimension with Dirichlet 
boundary conditions and with "Navier" ones; 

(II) Eigenfunction technique and nonlinear capacity method [61] are used in Section [3] 
to prove blow-up in finite time of solutions of these IBVPs with non-standard boundary 
conditions; 

(III) In Section m Henry's version of weighted Gronwall's inequalities are used prove 
global existence for the Cauchy problem in R N x R , for modified Kuramoto-Sivashinsky 
equations; 

(IV) In Section El we use two types of scaling blow-up techniques for global existence 
for the Cauchy problem and IBVPs. 

Our analysis embraces a number of 2m th-order iV- dimensional Kuramoto-Sivashinsky 
type models posed in R x R + . 




(1.2) 




for t > 
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In general, as customary in PDE theory, the following approach is classical: 



(1.3) local existence + a global a priori bound ==>- global existence. 



More precisely, in order to get (11. 3ft . we use the following intermediate blow-up step: 

(1.4) local existence + a global a priori bound no finite-time blow-up. 

Here scaling arguments to prevent blow-up are crucial. 

In this context, the problem of global solvability takes the following obvious negation 
form: 



1.5) global solvability = /9 finite-time blow-up of local smooth solutions. 



Of course, here we understand the a priori estimates in the corresponding function spaces. 
This proclamation is also relevant to the Navier-Stokes equations in IR 3 , which have a long 
mysterious history of uniquenes/non- uniqueness and blow-up singularity open problems. 

As rather unexpected but related application, in Section [6j we consider the Navier- 
Stokes equations 

(1.6) v t + (v • V)v = - Vp + Av, divv = in R N x R + , 

with bounded integrable divergence- free data v . It is worth noting that the convective 
term in the first equation in (11.61) for the velocity field has indeed a nonlinear dispersion 
mathematical nature as in the Kuramoto-Sivashinsky equation (II. ip . Therefore, it is 
rather natural to include the model (11.61) in the present context of semilinear KS-type 
equations. 

Blow-up self-similar singularities with finite energy do not exist. The idea that the clas- 
sical fundamental problem of unique solvability of (11.61) in IR 3 is associated with existence 
or nonexistence of certain blow-up singularities as t — > T~ , goes back to Th. von Karman; 
see [44]. Later on, in 1933-34, J. Leray [52], [54] proposed a mathematical question to 
look for blow-up in (I1.6P for N = 3 driven by the self-similar solutions of the standard 
dimensional type, with blow-up tim^3 T = 1, 

(1.7) v( x ,t) = ^= t w(y), p(x,t) = j±- t P(y), where y = 
Substituting (11.71) into (I1.6P yields for functions w and P a "stationary" system, 

(1.8) \w + \(y- V)w + (w- V)w = -VP + Aw, divw = in R N . 

During last fifteen years, a number of negative answers concerning existence of such 
non-trivial similarity patterns (11.71) . (11.81) in M 3 were obtained; see [131 E21 BO], and the 
most advanced and justifying negative answer in [41J. Let us note an existence result in 
[16] for N = 4. 

Nonexistence of Leray's similarity solutions (I1.7P and other local types of self-similar 
blow-up is a definite step towards better understanding of the singularity nature for 
the Navier-Stokes equations. Of course, this does not settle the problem of singularity 



1 Actually, Leray suggested also to look for a self-similar extension of the solutions beyond blow-up, 
i.e., for t > T using the same scaling similarity variables, [54] p. 245]. 



formation (or nonexistence of finite energy singularities which is more plausible), since 
there might be other ways for (11.61) to create singularities as t — > 1~ rather than the 
purely self-similar scenario ( 11.71) . This multiplicity question is discussed below. 

On a countable set of blow-up patterns with infinite energy. For infinite L 2 -energy, the 
blow-up (also called the enstrophy blow-up of vorticity) in the Navier-Stokes equations 
(11 . 6 p can occur even for N = 2. Such global blow-up described by von Kdrmdn solutions 
is explained in [321 Ch. 8], where a rigorous theory of such singularities was developed. 
Earlier history of such solutions can be found in pp. This blow-up creates a plane jet. It 
is worth noting that the blow-up behaviour is also not of self-similar form as t — > T~ , and 
it is given by a similarity solutions of a non-local first-order Hamilton- J acobi equation 
associated with such a flow. Moreover, and this is also crucial, that [3"2"l pp. 232-235] 

(1.9) there exists a countable set of such different blow-up patterns (N = 2). 

Similar solutions can be constructed for axisymmetric flows in cylindrical coordinates 
for the Navier-Stokes equations (11.61) in M 3 ; see an example in [H Ch. 7, § 3] and [63] . 
However, the mathematics of such blow-up patterns becomes more involved, and (11.91) for 
N = 3 demands difficult proofs [29] . 

Singular (blow-up) set has zero measure. There exists another classic direction of the sin- 



gularity theory for the Navier-Stokes equations that was originated by Leray himself [5] 
(details are available in [20J ) and in Caffarelli, Kohn, and Nirenberg [9]. It was shown that 
the one-dimensional Hausdorff measure of the singular (blow-up) points in a time-space 
cylinder is equal to zero. We refer to [551 159"] for further development and references. In 
particular, among other results including Leray's one in [51] , a refined criterion is obtained 
in [69], saying that, if t — 1 is the first singular (blow-up) moment for a solution v(x,t) 
of flTJD, then 

i 

(1.10) hm_j^—j;J J \v(x,t)\ 3 dxdt = +oo. 

t R3 

Performing the scaling as in ( 11.71) . 

(1.11) v(x,t) = -^=w(y,r), Z/=7=, r = -ln(l -t) -> +oo as t -> 1~, 
( 11.101) takes the same form 



+oo 

1.12) lfm e r J e ~ S (y |w(z,s)| 3 dy) ds = +oo. 



This means that, for the existence of a singular point t — 1, the solution of the rescaled 
equations, 

(1.13) w r + |w + i(y- V)w + (w- V)w = -VP + Aw, divw = in R N , 
must diverge (blow-up) as r — *• +oo in L 3 (R 3 ). 
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Thus, according to the criterion (j!.12p . t — 1 is not a singular (and hence regular) point, 
if the corresponding locally smooth solution of (11. 1 3 p does not blow-up as r — > oo in a 
suitable functional setting. Hence, the problem of global existence and uniqueness of a 
smooth solutions of the Navier-Stokes equations in M. 3 reduces to nonexistence of blow-up 
in infinite time for the rescaled system (I1.13p . In such a framework, this problem falls 
into the scope of standard blow- up/non-blow-up theory for nonlinear evolution PDEs. 
Countable sets of blow-up patterns in combustion problems. Actually, there are many ex- 
amples of countable sets of blow-up patterns for of much simpler reaction-diffusion equa- 
tions. Amongst them is the classic Frank-Kamenetskii equation (1938) [23] developed in 
combustion theory of solid fuels (also called solid fuel model), 

(1.14) u t = Au + e" in R N x R+ (N = 1, 2), 

for which there exists a countable set of different blow-up patterns. Rigorous mathematical 
theory of such blow-up patterns was known since the beginning of the 1990s (see e.g., 
[75| [76] ) and was developed by linearization in the inner blow-up region and nonlinear 
matching. A similar strategy to construct a countable set of blow-up patterns is applicable 
to the higher-order reaction-diffusion PDEs [26] 

(1.15) u t = -(-A) m u+\u\ p - l u (m>2, p > 1). 

where the analysis uses polynomial eigenfunctions and discrete spectrum of some related 
linear non self-adjoint differential operators [17] . 

On the other hand, the quasilinear ID counterpart of (11.141) with the p-Laplacian, 

(1.16) u t = {\u x \°u x ) x + e u (<7>0), 

is also known to admit a countable set of blow-up patterns [8] (see also a general discussion 
in [281 PP- 30-34]), but now, depending on a > 0, first few are self-similar, i.e., represent 
the case of nonlinear eigenfunctions (not linearized as above for (11.141) ). 
Evolution completeness as a necessary ingredient. It is a principal open problem to de- 
scribe the whole set of all possible blow-up patterns (if any) for the Navier-Stokes equa- 
tions (11.61) . Evidently, proving nonexistence of all the blowing up patterns, i.e., nonexis- 
tence of blow-up at all, will settle the fundamental problem of global smooth (and hence 
unique) continuation of sufficiently arbitrary solutions. Another natural possibility is to 
establish that, for given class of data, the orbits do not approach any of blow-up pattern 
scenario, so remain regular for all times. 

In this context, the problem of evolution completeness of the given countable set of 
patterns occur, meaning that these patterns exhaust all possible types of approaching 
the singularity for a fixed class of data. For linear problems, the evolution completeness 
follows from standard completeness and closure of eigenfunction subset of a linear operator 
or pencil in a fixed functional framework, so the evolution completeness does not have a 
separate meaning. For nonlinear problems, where "linear" notions of completeness and 
closure make no sense in general, the evolution completeness becomes key. 

There are some examples of evolution completeness of countable sets of linearized or 
nonlinear patterns (eigenfunctions) in parabolic asymptotic theory. For instance, the full 
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classification of blow-up sets for the Frank-Kamenetskii equation (j!.14p was performed 
by Velazquez [75] by actual proving the evolution completeness of the countable set of 
linearized blow-up patterns. Concerning nonlinear patterns, it seems that there exists a 
unique example of a proof of evolution completeness of such a countable set for the porous 
medium and p-Laplacian equations in ID or in radial geometry in M , 

u t = A(|w| m -V) and u t = V • (iVu^Vu), with m > 1; 

see [27], where the notion of evolution completeness was introduced. 
On the well-posed Burnett-type equations. Concerning our conclusion, as a straightfor- 
ward application of the scaling technique in (IV), we present a simple proof of the fact 
that a local smooth solution v(x,t) of (II. 6p . which is uniformly bounded in L P (M N ), with 

(1.17) p > N, i.e., p > 3 for N = 3, 

is uniformly bounded in L°°(M. ), so that such L p -solutions do not blow-up. The non- 
blow-up also is proved in the well-known critical case p = N = 2. The condition (11.171) 
is consistent with Leray-Prodi-Serrin-Ladyzhenskaya regularity L p,q criteria and other 
more recent results; see key references, history, details, and results in recent papers [SSI 
[201 ESI [25] . In addition, (11.171) re-derives Kato's class of unique mild solutions in M , [12] ; 
see details and key references in [25], [77] . 

This approach also extends to the related 2mth-order well-posed Burnett equations, 

(1.18) vt + (v ■ V)v = -Vp - (-A) m v, divv = in R N x M + , 

containing the higher-order diffusion operator — (— A) m v, with arbitrary m > 1. Then 
the regularity criterion similar to that in (11 . 17)) for m = 1 reads 

P > 2^1 fOT 771 = 1,2,3,.... 

On the other hand, we prove that for smooth fast decaying divergence-free L 2 -data v , 
finite-time blow-up is impossible in dimensions 

N < 2(2m - 1) (N = 2 for m = 1 is included, [531 [50]), 

so there exists a unique global classic bounded solution. 

2. Method of interpolation: global existence for the KSE 

2.1. A priori estimates. To demonstrate these classic approaches, we consider the one 
dimensional KSE in the following IB VP setting (for convenience, here D = ^-): 



(2.1) v t + D 4 v + D 2 v = \ Dv 2 , t>0, xe(-L,L), 

(2.2) v(x, 0) = vo(x), x e (-L, L), 
with either: 

(2.3) v = Dv = 0, x = —L, x = L, t > 0, or 

(2.4) v = D 2 v = 0, x = —L, x = L, t > 0. 
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A priori estimates. Multiplying (12.11) by v and integrating by parts over Q = (-L, L) 
with regard to (12.31) or (12.41) . we find 

(2.5) \f t j v 2 (x, t)dx + J \D 2 v\ 2 dx- J \Dv\ 2 dx = 0. 

n n n 

Due to ([ZSD or flU]), we have 

(2.6) J \Dv\ 2 dx < v 2 dx) 3 ( J \D 2 v\ 2 dx 



n 



E(t) = / v 2 (x,t)dx, 



Denoting by 



from (12.51) we get 

(J \D 2 v\ 2 dxfE^{t)- J \D 2 v\ 2 dx<\E{t). 
n n 

So 

(2.7) E(t)<E(t):=E(0)es. 

Next, integrating (12. 5p in t > 0, we see that 

t t 

,2,) J jm***«=J J\»?***-\m)+\E<!>). 

on on 
Due to (J2ISD or we have 

i t i 

0= D(vDv)dxdt= / / wL)\;da;dt+ / / |Dw| 2 dxdt 



on on on 

From this, it follows that 

-t r If 1 



(2.9) / / \Dv\ 2 dxdt<~ [ ! v 2 dxdt + \ [ [ \D 2 v\ 2 dxdt. 
Jo Jn ^ J Jq 2 J J n 

Thus, (12. 8j) implies that 

i t 

\j j\D 2 v\ 2 dxdt<\ J E{t)dt- l -E{t) + l ~E{$). 
on o 
Hence, by (12.71) one has 

t 

(2.10) J J \D 2 v\ 2 dxdt < 2£(0)e^ - E(0). 

o n 
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Then from (12.81) it follows that 



z 

(2.11) J J \Dv\ 2 dxdt < (2es - |)S(0). 



o n 



Estimate of J J v 2 \Dv\ 2 dx dt. Thanks to embedding inequality for C K we have 
o n 

\Ht)foo<CooJ \Dv\ 2 {x,t)dx 
h 

since = 0. Then 

2 



w 2 |£)w| 2 dx < Coof y |£to| 2 dz 
h 



Next 



|Di;| 2 dx) 2 < Jv 2 dx- J \D 2 v\ 2 dx = E(t) J \D 2 v\ 2 dx. 
n n n n 

Applying ( 12.1 Oft , we find 

t t 

on on 

(2-i2) } r 

< CooE(t) / / |D 2 w| 2 dxdr < Coo E{t)(2e^ - l)E{0) 
o n 

= Coo et(2e^ -l)£ 2 (0). 



Next, by using the above estimates and the linear theory of the parabolic equations of 
the 4th-order we get (here Qt = ^ x (0, T)) 

IK^II^ + lkll^ + ll^ll^) 

+\\D*v\\h (QT) + \\Dv\\h iQr) + HI HQt) 
{ ' ' <C(\\vDv\\ 2 LHQT) + \\v \\ 2 L2(n) ) 

<C 1 (e T E 2 (0)+e?E(0) + \\v \\ 2 L2m ) 
with a positive constant C\ independent of v and T > 0. 



2.2. Global existence. 



Theorem 2.1. For any v G L 2 (Q), both the IB VPs (TO- COp and (HOD, (TO) . Q 
nave solutions for any t > 0, satisfying inequality ( \2.13b . 

Proof. The ground a priori estimates (12. 7p . (12. 9p . (12.101) and (12.121) are obtained by 
multiplication of equation (12.11) by v. Thanks to this we use the Galerkin approach and 
these a priori estimates are valid for the Galerkin approximations {v m }. 

The passage to the limit v m — > v as m —>■ oo in the nonlinear term v(x, t)Dv(x, t) (x G 
Q C R, t > 0) is proved by standard arguments (see [5], [56]). As a result we get a weak 
solution. Then from (12. 12j) due to the linear theory of parabolic equations of the 4th order 
we obtain the final estimate (I2.13p . □ 

3. Method of eigenfunctions: blow-up of solutions 

We now show how, using a similar interpolation-eigenfunction technique, to derive 
sufficient conditions of blow-up for the KSE with special "boundary conditions" . 

3.1. Basic computations. Multiplying equation (12. lft by a function ip(x) belonging to 

C7 4 (R), we get 

(3.1) ± J vip dx + J v(D*ip + D 2 ip) dx = -± J v 2 D4>dx + B(v,tp), 

n n n 

where f2 = (0, L) C R and 

(3.2) B(v, ip) = [vD 3 ip - Dv ■ D 2 ip + (D 2 v + v)D^p - (D 3 v + Dv)^ + \ v 2 ^\ \. 
For the function 

tp(x) = 4>\(x) := \x — L\ x 

with A > 6, we have from (13. ip 

(3.3) 4 / ^ Ax = \ J v 2 ■ \x - L\ x - X dx - J v, h dx + *(.), where 

n n n 

IA-4 



(3.4) 

Denote by 



^ = A (A - l)|a; - L| A ~ 4 [(A - 2)(A -3) + \x- L\ 2 } and 
B (v) = -\ LV(0, t) + AL^ 1 ■ [(A - 1)(A - 2)L- 2 ]v(0, t) 
+L X [1 + A(A - l)L- 2 ]Dv(0, t) + XL x - 1 D 2 v(0, t) + L x D 3 v(0, t). 

J(t) = J v(x,t)\x- L\ x dx, tt=(0,L). 
n 
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Thanks to the inequality 

n a 

<\Jv 2 \x-L\ x - l dx + C x {L), with 
h 

C X (L) = A(A - 1)L^ [ (A - 2 ^- 3)2 + 2(^)(a-3) L 2 + _i_ L 4] ) 
from (13.31) it follows that 

dJ 



dt 

By using the Holder inequality 



<| J v^x-L^dx + Bo^-C^L). 



J 2 (t)<^ \ v 2 \x-L\ x ~ x dx, 



9. 

we get 

dJ 
dt 



(3.5) — > ^ j^J 2 (t) + B (v) - C X (L). 



3.2. General initial boundary value problem. We consider a non-standard IBVP for 
the KSE ( 12.11) that includes initial data 

(3.6) v(x, 0) = vo(x) 
and some general boundary conditions to be specified, 

(3.7) B Q (v) = h(t) at x = for t > 0. 

We understand a solution v of (12. ip . (13.61) . (13.71) as a weak solution in the sense of 
identity (13. 7\i with respect to any test function ip e C 4 (M) with an initial data in the 
sense 

v (x, t)ip(x) dx — > / v (x)ijj(x) dx as t — ► 0. 



Of course, the function v is assumed to belong to the function space such that identity 
(13. 7p makes sense. 

3.3. The blow-up results. As follows from (13.51) the blow-up phenomenon for IBVP 
flgZTJ) , fl3T6|) . fl3T7D depends on 

(3.8) H x {v,L):=B Q {v)-C x {L). 

It means it depends on the relationship between initial data vq(x) and boundary condi- 
tions. 
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Theorem 3.1. Let H\(v, L) > a 2 > for some A > 6 and L > 0. Then there is no global 
solution of §21) . flX6j) . 

Moreover, there is no global solution of (jff.ip . (15*. ?P /or x G (0,L) ; and 



> ^ tan(a/ct + Co), where 

k = (^±2±L-( 2+ ^)\ a = Va?>0, c = arctan^, J = J v (x)\x - L\ x dx, 

n 

and the blow-up time is estimates as: 



-co 



< - ■ 

00 (IK 

The proof follows immediately from (13.51) . 

Corollary 3.1. Suppose thatv(x,t) satisfies, 
v(0,t) = 0, 

Dv(0, t) + D 3 v(0, t) > const > Vt > 0, and 
Dv(0,t), D 2 v(0,t) > 0. 
Then there is Lq > such that all the assumptions of Theorem \2.1\ are fulfilled. 

Theorem 3.2. Let H\(v, L) > for some A > 6 and L > 0. Lei 

Jo := / v (x)\x - L\ x dx > 0. 



TTien there is no global solution of (\2.1\i . i\3.6} . I\3. 9j) 



Moreover, there is no global solution of (1ff.il) . ( Iff.gj) . Qff.7|) /or x G (0,L). In particular 
J satisfies 

J(t) > jztrv 
so that for the blow-up time T^ we have 

T ^ 1 

The proof follows from (13.51) . 

Theorem 3.3. Let H\(v, L) > —a 2 for some A > 6 and L > 0. Let 

Jo '■= f Vq(x)\x — L\ x dx > — . 



JTien i/iere is no global solution of (lff.il) . ( Iff.ffl) . (5.?) 



Moreover, there is no global solution of (Iff, ip . (\3.6} , i\3. ?p /or x G (0, L). In particular 
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Hence, the blow-up time satisfies: 

T 00 <^ln^^, where « = (^±^L-( 2 + A )) s onda=|a|. 
Again, the proof follows directly from (|3.5|) . 

Corollary 3.2. Lei fo(x) > c|x| M wif/i some // > and c > 0. Suppose that the boundary 
values v, Dv, D 2 v, and D 3 v at x = do not depend on t > 0. Then there exists A > 
such that the assumptions of Theorem \3.§ are fulfilled. 

4. Global existence and L°°-bounds by weighted Gronwall's inequalities 

4.1. A general KS-type model. In this section, we extend global existence approaches 
to more general 2mth-order parabolic equations of the KS-type. As it has been seen, the 
zero boundary conditions of the IBVPs under consideration reinforced the energy control 
and helped the global solvability of the KSE. It is then reasonable to consider the KS- 
type equations in unbounded domains suggesting "infinite propagation" with no spatial 
"obstacles" and bounds, which the Cauchy problem is most natural for. 

Thus, as a basic and typical model, we consider the Cauchy problem (CP) for the 
equation, which includes both stable and unstable linear diffusion terms as well as the 
"convection" : 

(4.1) v t = -(-A) 2 W (-A)WB!(M P ) in R N x R + [m = 21 = 2, 4, 6, ...), 
where p > 1 and Bi is the linear first-order differential operators in divergence form 
(4-2) B lU = I E (fc) d k D Xk u. 

For simplicity, the coefficients {d^} of the operator are assumed to be constant, though 
some (x, ^-dependence can be easily allowed. The classic KSE (12.1 p corresponds to 

(4.3) p = 2, 1 = 1, and N = 1 (d 1 = 1). 

We refer to (14. ip as to the modified Kuramoto-Sivashinsky equation (the mKSE). 
Thus, we consider for (14.1 1) the CP with bounded sufficiently smooth data 

(4.4) v(x,0)=v (x) in R N , v G L°°(R N ) n H 2m (R N ). 

For instance, for the main demonstrations of the techniques involved, we may always 
assume that vq(x) is smooth and has exponential decay at infinity. 

One of the key mathematical features of the mKSE (14.11) is that it admits multiplication 
by v in L 2 to get the a priori L 2 -bound. Namely, similar to related straightforward 
manipulations in Section [21 we have by the Holder inequality: 

, , Is IK*) II! = -l|AMt)ll! + / v(-A) l v dx < -\\^ l v(t)\\l + IKt)!| 2 ||Ayt)i| 2 

(4.5) J 

<ilK*)ll! =^ IK*)lll<INIiM for all t>0. 
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Another obvious higher-order model with a similar L 2 -control (14.51) is 

(4.6) v t = -(-A) m v + \v + B 1 (\v\ p ) in R N x R + (m = 2, 3, 4, ...), 

with the same convection-nonlinear dispersion and the linear unstable zero-order term 
+jV. On the other hand, \v can be replaced by —Av as in (ll.ip . However, for some 
2mth-order KS-type models, deriving L 2 -bounds can represent a hard problem, so we 
avoid using those in future application of the scaling and other techniques. 

Thus, as a principal issue of our analysis, the intention is to use this a priori estimate 
(14. 5 p to derive stronger uniform bounds on solutions v(x, t) in the sup x -norm for all t > 0. 
The developed scaling technique is rather general and applies to other higher-order KS- 
type models with a priori known L 2 , L p , Sobolev, or other types of weaker or stronger 
integral bounds on solutions. 

4.2. Fundamental solution. We will need the fundamental solution of the correspond- 
ing linear parabolic (poly-harmonic) equation 

(4.7) u t = —(—A) m u in R N xl + (m = 21), 
which has the self-similar form 

N 1 

(4.8) b(x,t) =t~^F(y), y = x/t^. 

The rescaled kernel F is the unique radial solution of the elliptic equation 

(4.9) BF = -(-A) m F + ^yVF + £F = in R N , J F = 1. 

F is oscillatory for m > 1, has exponential decay as \y\ — > oo, and satisfies, for some 
positive constants D and d depending on m and N [18], 

(4.10) \F(y)\< De~ d \y\ a in R N , where a = 6 (1, 2). 

4.3. Local existence and uniqueness of smooth solutions. As a standard practice 
(see [El [24] and [73j Ch. 15]), for sufficiently regular initial data, local in time existence 
and uniqueness of the classical solution of the CP ( 14.11) . ( 14.41) is studied (via Duhamel's 
principle) by using the integral equation, where we have integrated by parts once, 

t t 

(4.11) v{t) = bit) * v + J (-A) l b(t - s) * v(s) ds + J Bi b{t - s) * \v\ p (s) ds, 



where b(t) is the fundamental solution (14. 8p and is the adjoint operator, 

(4.12) B;(.) = -i£ ( *)A.*(*(0)- 

For more singular data Vq G L q , with q G (l,oo), the solution may not be classical and 
is then understood as a proper continuous curve u : [0, T] — > L q satisfying equation 
(14.111) . Such questions of local existence and uniqueness were first systematically studied 
by Weissler in the 1970s and 80s, [78], [79]. See the results in [791 PP- 87-90], which actually 
apply to 2m-th order equations like (14. ip . More recent results on local and global existence 
for higher-order parabolic equations such as (14.111) can be found in [21 [151 El I 
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In what follows, we are interested in global existence and behaviour of solutions, so we 
have assumed that data Vq are sufficiently regular, and then we will use the equivalent 
integral equation ( 14. lip for such purposes. 



4.4. Global existence by Henry's version of weighted Gronwall's inequality. 
Theorem 4.1. The Cauchy problem <\4-l[, fl^MI) has a global unique classical solution if 



(4.13) l<p<3 and p < p = 1 + 2(2 ^~ 1} . 

For the original KSE (12. ip with parameters (14. 3p . (14.131) is valid in dimension 

(4.14) N <6. 

Proof. Let us write (14. lip in greater detail, 

t 

N+21 



v(x,t)=t-*z j F(-)v (y)dy + J (t - s)~— ds j (-A) l F(-)v(y,s)dy 

(4.15) t 

+ J(t-sr^ds j BlF(^^)\vny,s)dy, (-) = ((^^r), 

RN 

where we mean in (—A) l F(z) and B*F(z) that the operators act in the ^-variable. Denote 

(4.16) V{t) = sup \v(x,t)\. 



Writing in ( 14. 15ft \v\ p = \v\ \v\ p 1 for p < 3 (p — 1 < 2), using Holder's inequality yields 

t 

\v(t)\ < sup \vq\ J \b(t)\ + J \A l b(t - s) * v(s)\ ds 



+ J \B\b{t- s)*{\v\ p ~ 1 \v\){s)\ds< Bup|u | ||F||i 



(4.17) 



t 

J 



C||A'F||i / V(s)(t-s)-™ds 



o 

t 



+ C\\BIF\\^_ j V(s)\\v(a)\\r\t-8)^ds, 



(4.18) where (3 = ^-2-tv(p-i) 



Am 
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Thus, by ( 14. 5 j) . we obtain weighted Gronwall's inequality (here j- — |) 

(4.19) V(t) <C + C J(t- s)-^V(s) ds + C J e^it - s)^V(s) ds, 



where, obviously, the last term on the right-hand side is key for boundedness of V(s). 
By Henry's estimates for such weighted inequalities [39, p. 188], it follows that V(s) is 
bounded for any t > 0, if (3 > 0, which is equivalent to the last inequality in (14. 13ft . □ 

4.5. On a double exponential L°°-growth by weighted Gronwall's inequality. 

As we have mentioned, for the KSE, it is principal also to establish the best estimate on 
the growth of global solutions for t 3> 1, to be compared with the exponential one ( 11. 21) . 

Let us see what kind of L°°-bound is guaranteed by the approach applied above. Con- 
sider the principal integral operator in ( 14.191) . where we skip all the constants C, 

t 

(4.20) V(t) < 1 + J e^it- s)^ l V{s)ds (/? > 0). 

o 

For p = 1, Henry's "discrete" proof in [391 P- 188] also gives the fact that solutions of 
such Gronwall's inequalities do not grow as t — > oo not faster than exponentially, which 
is fine according to (11.21) . 

(p—l)s . . 

With the presence of the multiplier e * in the kernel in (14.201) . this is not the case. 
Using the idea of Henry's proof, a certain estimate of the behaviour of V(t) for t ^> 1 
can be indeed obtained. Nevertheless, rather surprisingly, this will not be an exponential 
growth, which may naturally look most plausible. Namely, it is easy to see that an upper 
bound on the growth of solutions of (14.201) is doubly exponential in the sense that their 
solutions cannot grow as t —>■ +oo faster than the "supersolution" (corresponding to "=" 
in (TC2QD) 

(4.21) V(t) = exp { [ pfain + e] t^e 12 ^}, with arbitrarily small e > 0. 

In other words, the approach based on weighted Gronwall's inequalities, though proving 
global existence, supplies us with a rather non-realistic (in comparison with (11.21) and 
(I4.5P ) doubly exponential L°°-bound (I4.2ip on solutions v(x,t). We will then need in 
Section [5] to improve this bound via a different scaling approach. 

4.6. Application to a non-divergent equation. A similar technique being applied to 
the non-divergent diffusion-absorption equation 

(4.22) v t = -(-A) m v - Ivl^v in R N x R + , 

yields global boundedness of solutions in the subcritical Sobolev range [331 § 2], 
(4-23) Kp<p s = ^. 

Further refined applications of the scaling technique for (14.221) are given in [12] (note that 
Theorem 4.1 proved for p > ps therein in § 4 applies to sufficiently small solutions only). 
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5. Global existence and exponential L°°-bounds by scaling techniques 
5.1. Global existence in subcritical range: C^-scaling. Cf. Theorem I4.lt 

Theorem 5.1. The Cauchy problem (4-1), (4-4) has global unique classical solution if 

(5.1) Kp< Po = l + ^r± (m = 2Z). 

Thus, for the standard KSE (12.11) with parameters (14.31) . in comparison with (15. ip . we 
manage to skip the first assumption p < 3. 

Proof. Assume that the local classical solution v(x, t) of the CP blows up first time at 
some finite t = T. Then v(x,T~) is unbounded in L°°(M. N ), otherwise it can be extended 
as a bounded solution on some interval [T,T + e), with a sufficiently small e > 0, by using 
the integral equation (14.151) with a contractive operator in C(R N x [0,5]), 5 > small, 
equipped with the standard sup-norm. 

Thus, we argue by contradiction and assume that there exist sequences {t k } — > T~ , 
{xk} C M. N , and {Ck} such that 

(5.2) sup \v(x,t k )\ = \v(x k ,t k )\ = C k -> +oo. 

R N x[0,t k ] 

Using a modification of the rescaling technique in [33J, we perform the change 

(5.3) v k (x, t) = v(x k + x,t k + t) = C k w k (y, s), where x = a k y, t = a 2 k m s, 
where {a^} is such that the L 2 -norm is preserved after rescaling, i.e., 

(5.4) K(0)||2=|K(0)|| 2 a k = C k ^^0. 
Therefore, by ( 14.5ft . for all s for which w k (s) is defined, 

(5.5) \\wk(s)\\l = -4^2 v 2 k (x,t)dx <\\v \\ 2 2 e^ for all se[-l,^). 

k k J k k 

As usual, such a rescaling near blow-up time, in the limit k — > oo leads to the so-called 
ancient solutions (i.e., defined for all s < 0) in Hamilton's notation [38], which has been a 
typical technique of reaction-diffusion theory; see various form of its application in [671 [32]. 

Let, according to (15.11) . p < Pq. Then, substituting (I5.3P into equation (14.11) yields that 
w k satisfies (as usual, m = 21) 

(5.6) (w k ) s = -(-A) 2l w k + fx k (-A) l w k + z/fcBiKI^ in R N x R, where 

2rn , _ 2(2m-l) 

(5.7) n k = C k N -> 0, v k = al^OT 1 = C\ X N -> 0, k -> oo. 

We next perform backward shifting in time by fixing so > large enough (this is possible 
in the time-interval in (15.51) since a k — ► in (15.41) ). and setting w k (s) = w k (s — s ). Then 
by construction, we have that 

(5.8) \w k (s)\ < 1 and H'^fclh < C on (0, sq) 

are uniformly bounded classical solutions of the uniformly parabolic equation (14. TJ . By 
classic parabolic regularity theory [18, 24J, we have that the sequence {w k } is uniformly 
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bounded and equicontinuous on any compact subset from R N x (0, s ). Indeed, the nec- 
essary uniform gradient bound can be obtained from the integral equation for (15.61) . or by 
other usual regularity methods for uniformly parabolic equations. Note that compactness 
in any suitable weaker topologies (see e.g., [7H1 ES]) is also acceptable, since passing to 
the limit we arrive at a weak, and hence, classical solution of the limit (simpler) parabolic 
equation. 

Therefore, by the Ascoli-Arzela theorem, along a certain subsequence, W)~(s) — > w(s) 
uniformly on compact subsets from R x (0, sq). Passing to the limit in equation (15.61) and 
using that both scaling parameters satisfy fi k , v k — > 0, yields that w(s) is a bounded weak 
solution and, hence, a classical solution of the Cauchy problem for the linear parabolic 
equation (14.71) . 

(5.9) w s = —(—A) m w, with data |iZ>o| < 1 and ||iZ?o ||a — C. 

Using the Holder inequality in the convolution (see the first integral in (14.151) ) yields 

N 

(5.10) w(s ) = b(s ) * w \w(y,s )\ < (s )~^||F|| 2 ||it)o||a < 1, 

for all so 3> 1. Hence, the same holds for sup y \w k {y, sq)\ for k 3> 1, from whence comes 
the contradiction with the assumption sup y \v k (y, s )\ = 1. 

Thus, v(x,t) does not blow-up and remains bounded for all t > 0. □ 

The proof is entirely local, so the result holds for the Cauchy problem as well as for any 
other homogeneous basic IB VPs for (14. ip . where the boundary conditions cannot generate 
blow-up on the boundary themselves. 

5.2. Exponential L°°-bound. We next improve the double exponential L°°-bound in 

Theorem 5.2. The global solution of the Cauchy problem {4-1), (4-4) the subcritical 
parameter range (15.11) . p < p , satisfies 

(5.11) sup xeRN \v(x,t)\ < C eT°* forallt>0, where 7o = 2 nTpo-p) ( > i)' 
Proof. We assume that ( 15.21) holds for a monotone sequence {tk} —> +oo. By ( 14.51) . i.e., 

(5.12) \Ht k )\\ 2 < \\v Q \\ 2 e%, 

we perform scaling (15. 3p by taking into account the exponential factor in (I5.12p targeting 
a uniformly bounded rescaled solution in the sense of (15. 8p . This yields (cf. ( 15.41) ) 

(5.13) a k = Ct"e&, 

and eventually we arrive at the rescaled equation (15.61) . where 

Hi _™ mtu (2m— l)tt, 

(5.14) n k = a* = C/ew and v k = C p k ~ po e^r^ . 
Assume now that (15.111) does not valid so that, along the sequence {tk}, 

(5.15) Ck = Kke yotk , where K k — > oo as k — > oo. 
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One can see that then the rescaled equation (I5.6P contains the parameters (recall, 70 > i) 

_2_ 4-,,. _ i m 

(5.16) a k = n k N e" — tfc -> 0, fi k = a k 2 -> 0, and v h = <~ po -> 0. 

Hence, repeating the arguments of the proof of Theorem 15. 11 we obtain the limit problem 
( 15. 9ft . which does not support the assumed growing behaviour. □ 

5.3. Global existence for the critical exponent p = p : (T — i)-scaling. 

Theorem 5.3. The Cauchy problem fl^-ifl , d^-^D /ios a/so global unique classical solution 
in the critical case 

(5.17) p = Po = 1 + 2 J 2 I^l ( m = 2Z). 

Remark: on application of the C^-scaling. For p = p in (I5.6p . /i^ — > 0, but 
z/fc = 1, so that passing to the limit k — > oo, for the limit function iOfe(s) — > we 
obtain the KSE in x (-oo, 0) without the unstable diffusion-like term: 

(5.18) w s = -(-A) m w + BiM P , |w(s)| < 1, ||w(s)|| 2 < C, sup y 0)| = 1. 

Thus, we need to prove that such a solution defined for all s < is nonexistent. 

On the one hand, this looks rather reasonable, since in the class of uniformly bounded 
L 2 n iJ 2m -solutions, the PDE (15.181) in Mr x IR + is a smooth gradient dynamical system 
admitting positive definite Lyapunov function as in (14.51) . 

(5.19) i£lK s )H! = - / w(s)(-A) m w(s) = -\\D m w(s)\\l < 0. 

Hence, the only equilibrium is that is globally asymptotically stable in L 2 (M> N ), so that, 
in view of the interior parabolic regularity for bounded solutions, for any suitable initial 

data w Q G L°° n L 2 , 

(5.20) w(y, s) — > as s — >• +oo uniformly. 

On the other hand, this is not enough to complete the proof, since we need the convergence 
(15.201) to be uniform with respect to data satisfying conditions (15.81) . This will prove the 
actual nonexistence of a solution to ( I5.18p . but is not that straightforward. Therefore, 
for convenience, we choose another, but indeed related, scaling technique to prove non- 
blow-up for the critical exponent, which also emphasize other important aspects of this 
evolution PDE. 

Proof of Theorem \5.3[ Thus, assuming first L°°-blow-up at t = T~ , we perform the 
time-dependent (T — t)-scaling of the orbit {v(-,t), t G (0,T)} of (14.11) : 

(5.21) v{z,t) = (T-t)-°w(y,T), a = y = ^rt, r=-ln(T-t). 

Without loss of generality, we assume that x = is a blow-up point. Then w solves the 
following exponentially perturbed equation: 

w T = AO) + e~5(-A)V woeH 2m nL°°, 

(y5 ' 22 ^ A(w) = -(-A) m w- ±y- Vw - aw + B 1 \w\ p . 
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Note that the L 2 -invariance of the scaling (15. 16[) : 
(5-23) \\v(;t)\\ 2 =\\w(;r)\\ 2 . 
We next need some auxiliary results. 

Proposition 5.1. Assume that, along a subsequence {r^} — > +oo ; 

(5.24) w(y, Tfc) — > uniformly. 

Then t = T is not a blow-up time forv(x,t) (so the singularity is removable). 

Proof. Consider the sequence of solutions {wk(y, s) = w(y, r k + s)} with vanishing initial 
data in L°° according to (15.241) . We now use the good and well-developed spectral prop- 
erties [17] of the linearized operator B* in (I5.22p defined in the weighted space L%(M N ), 

where p*(y) = e~ a ' y ' /3 , (3 = , and a > is small enough: 

(5.25) B*=-(-A) TO -£y.V-aI, a(B*) = {A, = -a - ±, 1 = 0,1,2,...}, 

with eigenf unctions (generalized Hermite polynomials) that compose a complete and 
closed set. Therefore, according to classic asymptotic parabolic theory (see e.g., [57J), 
we conclude that for any sufficiently large k, 

(5.26) w k {y,s) ~ 0(e~ as ), s>l w(y, r) ~ 0(e~ aT ) , r > 1. 
Overall, taking into account scaling (15.211) . this yields: 

(5.27) v(x,t) ~ (T -ty a O(e- ar ) = 0(1) as t^T~, 
so that f(x, t) is uniformly bounded at t = T. □ 

Proposition 5.2. There must exist a subsequence {r^} — ► +oo such that 

(5.28) Ck = Tfc) 1 1 oo — *• °° as +oo. 
Proof. Assume for contradiction that 

(5.29) w(y,r) is uniformly bounded in M. N x ]R + . 

Then, similar to (I5.19p . we get for uniformly bounded and smooth solution w(y, r) that 

(5.30) iiHT)||3 = -||5™u>(r)||l + e-i||5 l «;(r)||l. 
Therefore, in order to avoid the L 2 -vanishing: 

(5.31) lk(r)|| 2 -^0 as r -> +oo, 

which by the interior parabolic regularity would imply (I5.24p as r — > +oo uniformly and 
hence no blow-up, one has to have that 

(5.32) f°° || D m w(r) || 2 dr < oo, i.e., ||D m w(r) ||| G #((1, oo)). 
Obviously, the convergence of the integral in (I5.32p implies that 

(5.33) D m w(-, r + s) ^ as r ^ oo weakly in Lf oc (R+; L 2 ). 
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Indeed, one can see that, by the Holder inequality, for any \ £ Cq°(WL n x R + ), 
| - JJ(-A) m w(r + s) X (s)\ = | ff(D™) 2 w(T + s) X (s)\ 

(5.34) 

= \JJ D m w{r + s)D m X {s)\ < {$j(D m w(T + s)) 2 y{jj(D m x(s)) 2 y ^0 

as t — > oo. Fixing a sequence {r k } — > oo and passing to the limit as r fc + s 
elude that, in view of the remaining interior parabolic regularity (recall that, regardless 
the "artificial degeneracy" (15.331) . equation (15.221) is uniformly parabolic), some equicon- 
tinuous subsequence {w{rk + s)} — > u)(s) uniformly, where w(y, s) is a smooth solution of 
the first-order Hamilton-Jacobi equation 

(5.35) w s = -^y -Ww - aw + B l \w\ p , w eL 2 nL°°. 

However, by classic theory of conservation laws [73] (and, actually, by the standard com- 
parison), all such solutions of (I5.35P are exponentially decaying: 

(5.36) w(y,s) =0(e~ as ) as s -> oo. 

This implies that there exists always a moment tq = r k + sq, with large enough k, r k , and 
So such that w(y,To) is arbitrarily uniformly. As in the proof of Proposition 15.11 overall, 
this implies no blow-up for v(x, t) at t = T. □ 

Finally, the result on global existence is completed by the following: 

Proposition 5.3. The Type if] blow-up solutions satisfying ( 15. 28b do not exist. 

Proof. We now apply in (15. 22ft the C^-scaling with {Ck} given in ( 15.281) . i.e., as in ( 15. 3ft . 
For p = p , this gives the following equation this gives the following equation for the 
rescaled sequence: 

(5.37) w k (y, t) = w(y k + a k z, r k + a 2 k m s) = C k w k (z, s), where w k (z, s) solve 

(5.38) w s = -(-A) m w + Bi|i&|p - a 2m (^ z ■ Vw + aw) + <e~ * * (-A) l w. 

Here, at s = 0, w G L 2 f] L°° and a k = C k 2/N . This is a uniformly parabolic equation 
with asymptotically small perturbations, so that on passage to the limit k — >• oo for this 
set {w k } of smooth solutions, one has to have in the limit that, along a subsequence, 

w k = w(y k + a k z, r k + a 2 k m s) — > W(z, s), where 

(5.39) 

W s = -(-A) m W + B 1 \W\ P in R N x R, W e L? n L°°, \\W(z, 0) ^ = 1. 

Note that s) ^ is an ancient solution, which is defined for all s < 0. At the same 
time, by construction, it is also a future solution, which must defined for all s > 0. Indeed, 
one can see that if W(s) blows up at some finite s = S~ > 0, this would contradict the 



2 The terms "Type I, II" were borrowed from Hamilton [38] , where Type II is also called slow blow-up. 
In reaction-diffusion theory [571(32], Type I blow-up is usually called of self-similar rate, while Type II is 
referred to as fast and non-self-similar. 
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Type II solution w(y, r) is globally defined for all r > 0. Thus, by scaling of the Type II 
blow-up orbit (15.281) . we arrive at the problem (I5.39p . which defines: 

(5.40) { heteroclinic solution W(z, s) ^ 0} = { ancient for s < 0} U { future for s > 0}. 

Let us more clearly specify the necessary properties of this mysterious and hypothetical 
heteroclinic solution W. Bearing in mind the possibility of multiplying (I5.39P by W and 
integrating by parts, which are again guaranteed by the convergence demand as in (I5.32p . 
we have the standard identity: 

+oo 

(5.41) lj- s \\W(s)\\l = -\\D m W(s)\\l<0, sGl =5> J \\D m W(s)\\ 2 ds < oo. 

— oo 

Hence, in the given class L 2 r\L°° classical solutions, (15.391) is a sufficiently smooth gradient 
system, and the only equilibrium that can be approached by such bounded orbits is W = 0. 
By passing to the limit s — > ±oo, we then obtain that uniformly 

(5.42) W(s) ^ as s -»• ±oo. 
This actually means that 

(5.43) the heteroclinic solutions W(z, s) ^ is a homoclinic of orbit . 
We next easily prove that 

Proposition 5.4. A nontrivial solution W ^ of (\5.39^i does not exist. 

Proof. It suffices to use the "ancient" part of the definition (15.401) . In view of the integral 
convergence in (I5.4ip at s = — oo, similar to (I5.24p we have that along any monotone 
sequence Sk — > — oo, we have that, for k 3> 1, W(sk + s) ~ V(s) uniformly on compact 
subsets in M> N x M, where V solves the Hamilton- Jacobi equation (a conservation law) 

(5.44) V S = B 1 \V\ P for s > 0, V = W(s k ) e L 2 n L°° . 

Since in the class of smooth solutions, V(s) decays fast (see e.g., [701 E3])> we have that 

(5.45) ll^ / (0)|| oo > ||y(s)||oo for 

so that the same is true for W(sk + s) provided that k ^> 1. In fact, (15.451) means 
that using proper theory of conservation laws (15.441) implies that a nontrivial L 2 H L°° 
ancient solution does not exists for such gradient system, since by the convergence in 
(15.411) . equation (I5.39P does not have a sufficient mechanism for growth of solutions from 
||W(— 00)1(00 = to || VK(0) ||oo = 1. To justify this more clearly, given an ancient solution 
W(z,s) of (I5.39p . we perform the standard scaling: 

(5.46) w^,s) = A-^(^,f), A>0, 

where W\ is also an ancient solution of (I5.39P for any A > 0. We now pass to the limit 
A — > + by using these simple properties: 

\\W x (0)\\oo = \- a ^ +00, ||W A (s)||l = \\W(f)\\ 2 -+ c > 0, 



(5.47) 

and $f s \\W x {s)\\r 
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and § £ ||W A ( S )||| = _A||D m H/(f)||i -> 0. 



Therefore, in the limit A — > 0, we have to have a nontrivial L 2 -solution V(z, s) of (I5.44p . 
which blows up as s — > _ . Obviously, this contradicts the Maximum Principle for this 
first-order conservation law, [73l Ch. 16]. 

In other words, we have that, for the dynamical system in (15.391) . 

(5.48) the unstable manifold of the origin is empty. □ 
This also completes the proof of Theorem 15.31 □ 

5.4. On uniform bounds in other KS-type models. For the non-divergent equation 
(14.221) . a similar scaling techniques yields that [33j, in the Sobolev range (I4.23p . solutions 
are uniformly bounded, i.e., 

(5.49) \v(x,t)\<C in R N x R + . 

For the divergent mKS-type equations without the unstable backward diffusion term, 

(5.50) v t = -(-A) m v + B l {\v\ p ) in R N x R + (m > 2), 
we easily extend the result as follows: 

Theorem 5.4. The Cauchy problem (\5.5Uty . i \4-4ty the parameter range (\5.1\i has a 
global unique classical solution, which is uniformly bounded, i.e., ( 5.4fy holds. 

Proof. Again, it suffices to consider the case when (15. 2p holds for some sequence {tk} — > 
+oo. Then the same proof leads to the contradiction with the hypothesis that {C&} — > 
+oo, and hence v(x,t) is global and uniformly bounded. □ 

5.5. On a generalization: higher-order nonlinear dispersion. In order to extend 
application of our final approach, we next briefly discuss mKS-type equations with a 
third- order (also odd) nonlinear perturbation (dispersion) of the form 

(5.51) v t = -(-A) m v - ABx(H P ) in R N x R + (m > 2). 
Writing this PDE in a pseudo-parabolic form, 

(5.52) Pv t = (-A)™-^ + B 1 (\v\ p ), where P = (-A)" 1 > 0, 

and multiplying by v in L 2 (R N ), we observe that, instead of a uniform L 2 -bound, we are 
given an a priori if ~ Abound: for uniformly bounded data 

(5.53) v E L°°(M. N ) n H 2m (R N ), 
the following holds: 

(5.54) ||v(t)||_i < ||v ||_i for t > 0. 

Here, for simplicity, we assume that vq(x) has also exponential decay at infinity, so v(x, t) 
does for t > 0. In (1 5 . 5 2 [) and later on, we define (— A) _1 w = g in a standard manner: 

(5.55) Ag = —w in K , g(x) — >• as x — > oo. 
For the solvability of this problem we shall always assumes that 
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(5.56) J wdx = 0. 

R N 

Clearly this property holds for the divergent equation (15.511) with exponentially decaying 
solutions. 

Theorem 5.5. A unique solution v(-,t) 6 iJ _1 (lR Ar ) for t > of the Cauchy problem 
(\5.51\) . i \4-4h ([S uniformly bounded, i.e., ( |5.^ffi ) holds, in the parameter range 



(5.57) Kp<p„ = l + M 

Proof. Local existence of classic solutions for (I5.5ip also follows from the equivalent 
integral equation such as (14. lip , with no second unstable term by replacing 

B* ^ -B*A. 

One can see that we still obtain a locally integrable in t > kernel that defines the 
operator being a contraction in C(R x [0,5]), with sup-norm. 

To prove global and uniform boundedness, we use the same scaling as in the proof of 
Theorem 15. 1[ where, instead of (15.41) . keeping the iJ _1 -norm of w k yields 



2 

' N+2 



(5.58) H(*k)ll-i = K(0)||_i =► a k = C k 

Then, for p < p , we arrive at the corresponding equation such as (15.61) . where 



-. 2(2n 

,2m-3/op-l _ " 



(5.59) v k = a* n - A C{r' = Cl N+2 ^0 (ji k = 0) 

as k — > oo. By passing to the limit k — ► oo, we again arrive at the purely poly-harmonic 
flow as in (15.91) that cannot support the necessary properties of the sequence. As usual, 
this analysis applies twice: (i) {t k } — > T~ < oo, to prove that no blow-up occurs, and (ii) 
{tk} ~^ +oo, to prove that the solution is uniformly bounded. □ 

The critical case p = po can be also studied in similar lines as above. Note in addition 
that is the unique globally asymptotically stable equilibrium of the smooth gradient 
system (15.511) in the corresponding class of regular solutions that are uniformly bounded 
in L°° n H~ x . 

Further generalizations including odd higher-order nonlinear dispersion operators are 
straightforward. 

5.6. On blow-up in divergent models. Phenomena of finite-time blow-up in such 
semilinear models with divergent operators is most well-known for unstable limit Cahn- 
Hilliard equation 

(5.60) u t = -A 2 u- A(H p -V) in R N x R + (p > 1). 
Blow-up solutions have the standard self-similar form 

(5.61) u(x,t) = (T-ty^f(y), y = x/(T-t)*, 
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where / solves the elliptic equation 

(5.62) _A 2 /-A(|/rV)-!y-V/-_i_/ = in R * 

This equation admits a complicated set of solutions. For instance, for N = 1 and p = 3, 
it has a countable set of different profiles / that describe various types blow-up patterns 
(15.611) . See [21] for earlier related references and further results. The case p = 2 was 
earlier considered in [I]. Concerning Leray's scenario of similarity blow-up in (15.601) in 
both limits t — > T T , see [29] and references therein. 

We are not aware of reliable traces of standard blow-up for KS-type PDEs such as (14.11) 
or (15.511) with odd-order nonlinear dispersion terms in the supercritical range p > po. 
Therefore, we do not know whether conditions (15.11) and (I5.57P of global solvability reflect 
the actual evolution properties of the PDEs under consideration, or are sometimes purely 
technical. In other words, one then faces another fundamental problem on construction 
of blow-up patterns for L 2 -bounded solutions. Regardless a sufficiently strong progress in 
understanding of formation of blow-up singularities in various nonlinear PDEs achieved in 
the last twenty five years (see references and results in the monographs [28l [32l [6T| [67 ] fTT]) . 
for higher-order equations, there are still several fundamental open problems in identifying 
admitted structures of blow-up patterns. 

6. On L°°-bounds for the Navier-Stokes equations in R n and well-posed 

Burnett equations 

6.1. A classical fluid model in R . Consider the Navier-Stokes equations (11.61) with 
given bounded L p -data v . In order to apply our scaling argument, we use the fact that a 
classical bounded solution v(x, t) can be locally extended by using the integral equation 
that is similar to (I4.15p . Existence of such a local semigroup of smooth bounded solutions 
for the NSEs is well known; see Majda-Bertozzi [58] . 

Let us present some comments that will be useful for the Burnett equations f II . 1 8 j) with 
m > 2. Taking the fundamental solution (14. 8p for m = 1 with the rescaled Gaussian 



(6.1) F(y) = (4tt)-t e~ 

we consider (II. 6ft with h = — Vp cLS cL system 



(6.2) 



t t 
v(t) = b(t) * v - J b(t-s)* [(v • V)v](s) ds + J b(t - s) * h(s) ds, 



t t 
div b(t) * v - J div b(t - s) * [(v • V)v](s) ds + J div6(t - s) * h(s) ds = 0. 



As usual, the second equation in (16.21) is the one for the pressure corresponding to the 
solenoidal vector field v. Observe that, due to the exponential decay of the Gaussian 
(16. ip . the first equation contains the operator in v that contractive in a bounded closed 

24 



subset M 5 of C([0, 5], C\R N )), where 5 > sufficiently small, with the sup-norm. Indeed, 
assuming that in Ms, 

|v| < C and |Dv| < C, 

we can use the possibility of differentiating in x the equation once to control D\. For 
fixed vectors v and h(s), the contractivity of the principal "convective" operator 

t 

(6.3) JV(v) = fb(t-s)* [(v V)v](s)ds 

o 

in for small 5 > is then straightforward. Note that the standard Picard iteration 
scheme for (16. 2p can be put into the probability framework by using the linear semigroup 
associated with 3D Brownian motion; see a survey [77] for these and other details related 
to the Navier-Stokes equations. We stop at this moment discussing the local interior 
regularity theory for the Navier-Stokes equations and refer to [TQl E7J [58] as a guide for 
detailed developments in this direction. 

As a standard classic alternative way for local applications, the pressure is excluded 
from the NSEs 

(6.4) vt = H(v) = Av - P (v • V)v, where Pu = u- VA _1 (V • u) 

is the Leray-Hopf projector onto the solenoidal vector field. Using the fundamental solu- 
tion of A in Mr, N > 3 (cr^ is the surface area of the unit ball B\ C M ) 

(6.5) b N {y) = - (N J- 2)aN j^pn, where a N = f^, 

the operator in (16.41) is written in the form of Leray's formulation [58j p. 32] 



H (v) = Av - (v ■ V)v + C 3 / rf^, tr(Vv(z, r)) 2 dz, 
(6-6) I 

where tx(Vv(z, r)) 2 = v % v i and C N = ^ > 0. 

Equivalently, the nonlocal parabolic equation (16 .4p is known to induce a local semigroup 
of smooth solutions. 

As usual, our intention is to show that local sufficiently smooth solutions cannot blow-up 
under a certain L p -type constraint. We next clarify the conditions that prevent finite-time 
blow-up of solutions in L°° and exclude Leray' similarity (ll.7p or any other. 

It is a classic matter that the L 2 -norm is natural for (11.61) . since after multiplication by 
v, the convective and pressure terms vanish on sufficiently smooth functions v(x,t) with 
fast decay at infinity, 

(6.7) ((v-V)v,v) = and - (Vp, v) = (p, V • v) = 0. 
Therefore, on smooth solutions, 

t 

(6.8) ||||v(t)|| 2 = -||L»v(t)|| 2 =► ||v(t)|| 2 + 2|p V (t)|| 2 dt<|| V o|| 2 , t>0. 

o 
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Actually, the estimate in (16.81) is the energy inequality for Leray-Hopf weak solutions of 
fll.6P ; see e.g., [HS] and references therein. 

Nevertheless, an L 2 -bound is not sufficient to control non-blowing up property of so- 
lutions, and this is the origin of extensive mathematical research in the last fifty years. 
Recall that global weak solutions of (11.61) satisfying 

u E L°°(R + ; L 2 (M 3 )) n L 2 (R + ; H 1 ^ 3 )) 

were constructed by Leray [54|, and Hopf [40j in 1951. 

6.2. Application of blow-up C^-scaling in the supercritical case p > N. Namely, 
as a first simple constraint, which is inherited from our previous study of the KS-type 
equations, we assume that 

(6.9) ||v(t)|| p <C for t>0 (p>2). 

Then, seeking L°°-bound and hence assuming (15.21) . we perform the scaling (15. 3p . where 
we impose the preservation of the L p -norm of the rescaled function, i.e., 

p_ 

(6.10) Hvfcllp = ||w fe ||p =*> a k = C k N -> 0, and h k = a\. 

As usual, we next perform passage to the limit in the NSEs. This can be done in the 
framework of the original model (11.61) . as well as of the nonlocal parabolic representation 
(16. 4p . which we actually do. Note that passage to the limit in the integral term causes no 
difficulty for our sequences of uniformly bounded smooth solutions {w^}, where Ascoli- 
Arzela classic theorem [4lH Ch. 2] applies. 

We then obtain the following rescaled equations for w = w k (y, s): 

(6.11) w s + u k P(w • V)w = Aw, where v k = C* - * 7 ^ for p > N. 
Next, after time shifting, s i— »■ s — s , the solutions and data satisfy (cf. (15.81) ) 

(6.12) \w k \ < 1 and ||wfe|| P < C. 

By regularity for uniformly bounded sequence {w fc (s)}, we conclude that there exists its 
partial limit w(s) satisfying the solenoidal heat equation 

(6.13) w s = Aw. 

Of course, this is an equivalent form of writing the nonstationary linear "convectionless" 
Stokes system, 

(6.14) w s = -Vq + Aw, div w = in R N x R + , 

with data as in (I6.12p . It is an exercise to check, by using the integral equation with the 
kernel as in (I4.15p . that this problem for (I6.13P or (I6.14p admits a solution that exhibits 
a power decay for s ^> 1. Namely, an estimate similar to (15.101) can be obtained with a 
different exponent of the decay rate. 

In other words, under the assumption (16.91) . the Navier-Stokes system does not have a 
mechanism to create any L°° blow-up singularities. We thus arrive at the following: 
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Proposition 6.1. Under the given assumptions, LP -solutions, with p > N > 2 in ( 16.91) of 

the Navier-Stokes equations ( li.ffi) do not blow-up, and, moreover, are uniformly bounded: 

(6.15) \v(x,t)\<C in R N xR + . 

The last bound (16.151) is proved as in Theorem 15.21 by assuming that {t k } — > +oo. 
One can change the functional space in estimate (16. 9p to get L°°-bound. 

6.3. Critical case p = N = 2: an example of application of the (T — t)-scaling. 

In the critical case p = N, we have that = 1 in (16.111) . so that in the limit k — > +oo 
we arrive at the same Navier-Stokes equations, but now with uniformly bounded data and 
solutions in both L 3 (R N ) and L 00 (M Ar ). Moreover, by passing to the limit Sq = s^j — > +oo, 
we actually deal with the following class of solutions: 

(6.16) w(s) G L°° fl L N for all s > 0: |w(s)| < 1, ||w(s)||jv < C. 

In other words, using the scaling in the critical case p = N, we eventually get into the spe- 
cial class of solutions (16.161) . so a key restriction of the Navier-Stokes equations is achieved. 
Obviously, no blow-up and other singularities are available in the class (16.161) . Thus, in 
the critical case, in the class (16. lb]) , the Navier-Stokes equations induce a smooth gradient 
dynamical system with the positive Lyapunov function as in (16.81) which is strictly mono- 
tone on such nontrivial solutions. Hence, this admits the unique globally asymptotically 
stable trivial equilibrium 0. 

However, for p = N, proving non blow-up of solutions leads to a hard problem of nonex- 
istence of suitable ancient solutions. Similar to the previous KS problem, we demonstrate 
an example of application of the (T — i)-scaling to get the result in the critical case 
p = N = 2 of the obvious particular interest. 

We begin with Leray's blow-up scaling [M] for (16. 4p by setting 

(6.17) v(x,t) = 7 fow(j/,r), y= 7 §= t , r=-ln(T-t), 
to get the following rescaled equation: 

(6.18) w T + P(w ■ V)w = B*w, where B* = A-§2/-V-§I. 
Here, B* is the adjoint Hermite operator with the discrete spectrum 

(6.19) a(B*) = {X k = -\ - f , k = 0,1,2,...} 

and a complete and closed set of eigenfunctions being finite solenoidal Hermite polyno- 
mials; see details in [3H Append. A] and [30l § 2] , where all the aspects of the functional 
settings of B* in the weighted space L 2 pt , p*(y) = / 4 , can be found. Note that scaling 
(16.171) implies the following behaviour of the rescaled solution: 

(6.20) || w (r)||^e^^| iv=2 = 0(l) as r - +oo. 

We next follow the same scheme as in Section [5j 

(i) First, the L 2 -conservation such as f!5.23j) holds for w(y,r); see ( I6.20p . 

(ii) Second, assuming (15.241) yields by using the spectral properties in ( 16.191) and scaling 
(BOH (cf. (1537]) ) the non-blow-up in L°° at t = T of v(x,t). 
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(iii) Thirdly, we arrive at (15.281) for w(y, r), since its uniform boundedness will lead to 
L°°-decay by the identity 

(6.21) ^l|w(r)||H-/|Vw| 2 <0. 

This V-property of the dynamical system means that, in the class of bounded smooth 
solutions, the only equilibrium is trivial w = 0, so that the uniform stabilization to it 
again guarantees non-blow-up. 

(iv) Fourth, assuming (I5.28p . we C^-scale such an orbit to get a convergence to a nontrivial 
0-homoclinic solution as in (15.391) : 

oo 

(6.22) W s + P(W • V)W = AW in R N x E, / ||VW(s)|||ds < oo, 

— oo 

where, thus, W(s) G L 2 n L°°, and, by fl6T2T|) . W(s) E H 1 for all s 6 R. Otherwise, 
if W H 1 , then (16.211) will guarantee fast decay of the L 2 -norm of w(-,r) and hence 
v (-,t) as t — > T~ , which is contradictory. Finally, using the converging integral in (16.221) 
and passing to the limit as + s — > — oo leads to smooth small solutions of the Euler 
equations (EEs): 

(6.23) V S + P(V- V)V = 0, V(s)GL 2 nL°°. 

Then the void conclusion similar to (15.481) remains valid, provided that smooth L 2 - 
solutions of the EEs (16.231) decay in L°° sufficiently fast, which is true since this is a 
smooth gradient system (recall that we are obliged to treat the simpler critical case 
p = N = 2 only); see [73J Ch. 17] and surveys [31 [H] for further details. This shows a 
potential correct way to treat the relations between singularities in the NSEs and EEs, 
where the absence of some of those for the latter ones implies the same for the former. 

Of course, our analysis of global existence of classical bounded solutions in the critical 
case p = N = 2 just reflects the classic Leray's (1933, the CP) [53J and Ladyzhenskaya 
(1958, IBVPs) [1H H9] (see also [50]) existence-uniqueness results for N = 2. In this 
connection, it is worth mentioning another new proof for N = 2 in Mattingly-Sinai |59j . 
which is based on using advanced Fourier Transform techniques^. 

6.4. 2mth-order well-posed Burnett equations. As a natural extension, similar to 
KS-type problems in Sections H] and we consider the 2mth-order model (jl.lBp . Ac- 
tually, for m > 2, such models are not that formal and have been known to appear as 
Burnett's equations on the basis of Grad's method in Chapman-Enskog expansions for 
hydrodynamics. Unlike (I1.18p . the original Burnett equations are ill-posed, as backward 
higher-order parabolic equations having a wrong sign at the diffusivity operators. Namely, 
Grad's method applied to kinetic equations yields, in addition to the classic operators of 
the Euler equations, other viscosity parts, as follows: 

oo 

v t + (v • V)v = £ £ 2jl+1 A"(/i ri Av) + ... = e(/i Av + £ 2 // 1 A 2 v + ...) + ... , 

n=0 



3 Actually, the proof therein looks not that "elementary" as the title of [55] suggests; we pretend that 
our approach is more elementary, though this is not a proper point for arguing. 
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where e > is essentially the Knudsen number Kn; see details in Rosenau's regularization 
approach, [BB]. In a full model, truncating such series at n = leads to the Navier-Stokes 
equations ( 11. 6ft (with fiQ > 0), while n = 1 is associated with the Burnett equations. 
These are ill-posed since, by expansion, \i\ > 0, so a backward bi-harmonic flow occurs, 
etc. We will refer to (I1.18P for m > 2 as to the well-posed Burnett equations. 

Note also that Burnett-type equations with a small parameter appeared as higher- 
order viscosity approximations of the Navier-Stokes equations is an effective tool for 
proving existence of their weak ( "turbulent" in Leray's sense) solutions; see Lions' classic 
monograph [561 § 6, Ch. 1]. 

For the system (11.181) . we use the following parameters of scaling: 

_p_ 1 p(2m-l) 

(6.24) a k = C k \ b k = a 2 k m , and D k = C k + N . 

This gives the parameter of the convection term in the analogy of 06. lip as: 

(6.25) v k = & k 



p(2m-l) 



N 
k 



Hence v k — > as k —>■ oo under the following hypothesis: 

Proposition 6.2. Under the given assumptions, LP -solutions of the well-posed Burnett 
equations <\ 1. 18§ . with 



(6.26) p > Po 



N 



2m- 1 



in ( I6'.^p . do not blow-up, and are uniformly bounded, i.e., (16.151) holds. 

The scaling analysis can be applies also directly to the locally smooth integral nonlocal 
parabolic flow similar to (16.41) . The critical case p = po is then treated by an additional 
use of the (T — t)-scaling, where some technical difficulties may occur. 

6.5. Well-posed Burnett equations: no blow-up for N < 2(2m — 1). This is a 
simple consequence of the previous scaling analysis. We recall that, for (11.181) . L 2 -norm 
of v(i) does not increase with time, so, for smooth solutions, the analogy of (16. 8p holds. 

Proposition 6.3. Let v G L 2 (R N ) r\LL 2m (M, N ) be divergence free. Then there exists the 
unique global bounded smooth solution of the well-posed Burnett equations (li.igp if 

(6.27) iV<2(2m-l). 

Indeed, substituting p = 2 into (16.251) yields that 

(i) v k — > in the subcritical case N < 2(2m — 1), where the proof is completed in similar 
lines and causes no extra difficulties. 

(ii) In addition, v k = 1 in the critical case 



(6.28) 



N = 2(2m - 1), 

29 



where the dynamical system is a gradient one in the class of smooth bounded solutions. 
However, the proof of global existence according to the corresponding scaling quite similar 
to that in (IBTTTIh 

(6.29) v(x,t) = (T-ty 2 -^w(y,T), y = r = -ln(T-t), 
will lead to the following rescaled equation: 

(6.30) w T + P(w V)w = B*w, B* = -{-A) m - ^y ■ V - ^ I . 

Therefore, a proper spectral theory of generalized solenoidal Hermite polynomials as eigen- 
functions of B*, with the spectrum 

(6-31) <r(B*) = {\ k = -2^1 - A fc = 0,1,2,...}, 

is necessary. This is developed along the lines presented in [17], where the scalar 2mth- 
order case was under scrutiny. Eventually, this will allow to treat (we do not guarantee 
that all steps are going to be simple as in Section 16.31) the global existence of bounded 
classical solution in the critical case ( 16.28ft . As we have seen, for m = 1, (16.281) reads 
N = 2 and reflects classic Leray's [53] and Ladyzhenskaya's [HJ HHl [50] results. 

As is well-known, existence or nonexistence of L 2 -bounded L°°-blow-up patterns in the 
complement to (16.271) . (I6.28f) range 

(6.32) N > 2(2m - 1), or JV > 3 for m = 1, including the crucial 3D N = 3, 

comprises the core of the Millennium Prize Problem for the Clay Institute; see Fefferman 
[22] . Concerning possible structures of blow-up patterns in the NSEs, see discussion and 
review of a large amount of existing literature in a [30] that seems to be a most recent 
survey in the area. 
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